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Abstract 

We make a complete analysis of radiative symmetry breaking in the MSSM and its SU (5) 
extensions including low- and high-energy threshold effects in the framework of the two-loop 
renormalization group. In particular, we consider minimal SU (5) , the missing- doublet SU (5) , 
a Peccei-Quinn invariant version of SU(5) as well as a version with light adjoint remnants. 
We derive permitted ranges for the parameters of these models in relation to predicted a s 
and Mq values within the present experimental accuracy. The parameter regions allowed 
under the constraints of radiative symmetry breaking, perturbativity and proton stability, 
include the experimentally designated domain for a s . In the case of the minimal SU(5), 
the values of a s obtained are somewhat large in comparison with the experimental average. 
The missing- doublet SU(5), generally, predicts smaller values of a s . In both versions of the 
missing- doublet, the high energy threshold effects on a s operate in the opposite direction than 
in the case of the minimal model, leading to small values. In the case of the Peccei-Quinn 
version however the presence of an extra intermediate scale allows to achieve an excellent 
agreement with the experimental a s values. Finally, the last considered version, with light 
remnants, exhibits unification of couplings at string scale at the expense however of rather 
large ot s values. 
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1. Introduction 



Supersymmetric unification, in the framework of supersymmetric grand unified theories 
(SUSY GUT's) [0, or Superstrings 0, is in good agreement || f§] with the low energy values 
of the three gauge couplings, known to the present experimental accuracy, as well as with 
available bounds on the stability of the proton. The effective low energy theory resulting from 
such a framework is a supersymmetric SU(3) C x SU(2)l x U(1)y, model with softly broken 
supersymmetry. The simplest model of that class is the minimal extension of the standard 
model (MSSM) |TJ. A most appealing feature realized in the MSSM is the breaking of the 
electroweak symmetry through radiative corrections |J. The Higgs boson mass squared 
parameters, although positive definite at high energies, are radiatively corrected, as can be 
most easily studied by the use of the renormalization group, yielding a negative mass squared 
eigenvalue at low energies which triggers electroweak symmetry breaking J7|. 

In the present article we study the radiative breaking of electroweak symmetry in the 
framework of SU(5) SUSY GUTs. In addition to the standard low energy inputs («bm, 
Gp, M z , ...) and the soft breaking parameters (M , M 1 / 2 , A ), we have the thresholds 
of the superheavy particles parametrized in terms of at least two more parameters (Ms, 
M Hc , ... ) . Our output includes the strong coupling a s and the unification scale M G , as well 
as the complete spectrum of new particles. The predicted strong coupling values can depend 
strongly on the high energy thresholds. Thus, our analysis discriminates between the various 
GUT models. 



Our basic low energy inputs are the boundary values of the gauge couplings in the DR 
scheme |[| a,\ = ai(M z )\jy^ and d 2 = a.2(M z )\-jj^. Their values can be determined in terms 
of the Fermi constant G F = 1.16639 x 10 5 GeV" 2 , the Z-boson mass M z = 91.1884 ±0.0022 
GeV ||, the electromagnetic coupling a^ M = 137.036, the bottom quark mass mj = 5 GeV, 
the tau mass m T = 1.777 GeV, and the top quark mass. We can write down the following 
formulas, 

«i = g«EA/Cos 0(1 - A 7 + ln(— )) (1) 

*-l -1 • 2/j/i a i a EM , ( M S,, / n x 

« 2 = a EM sm 6{l- A 7 + ^— ln(— )) (2) 



where A 7 =0.0682±0.0007 |10[ includes the light quark and lepton contributions, and the 
DR mixing angle is given by 

sin6=-\l— 1 j= 2^ (3) 

2 1 1 V2G F M Z {1- Ar) J 1 V ; 

The quantity Ar will be given below. The scale Ms appearing in (1) and (2) is not a physical 
scale but a convenient parametrization for the contribution of all sparticle and heavy particles 
(W,t } H+), defined as 

MI 448 8 1 1 2 2 111 1 11 

Mw 7 M t 9 Ml+M? M? M* n Mge M? Ml Mf M\ M| M? M~. „ Ml Ml Ml 

W t H+ ti t 2 "1,2 u 1 2 J, ft 2 rfi 2 d? o T1 T2 ei > 2 e i,2 XI X2 

M 5 = M : ^ (4) 

Mi 
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An additional useful parametrization scale Ms, relevant to the strong coupling constant a s , 
is 



M?MfM?Ml ML Mi Mi Ml M§ c Mi 
M s = it ■ '■ (5) 



Mi 



The quantity Ar appearing in Eq. (3) can be written as [JTT 



Ar = A 7 - <f± log Ms _ n^S) + IWO) + A 5M + ^ + (6) 
27T M z Mf A% 



FI^z and Ilvyiy stand for the Z and W self-energies calculated using dimensional reduction. 



The quantity Asm stands for Standard Model vertex+box corrections and is given by [12 



A ol em . log cos 2 9 w 7 5sm 2 6»vy . 2 , 3cos 2 9 w 

= i™^ {6 + — 2 — sm e{5 ~ -2^f )]} (7) 

In the last formula, by definition cos 2 6w = -jp~- The pole mass of W-gauge boson in (4) is 
related to Mz by = Mf pcos 2 9 where the p parameter is given by, 

p = 1 WW fai W ^ ^ Z ^ + 2 — loop finite corrections. (8) 

M w M z 

We have explicitly written in (6) the 2-loop corrections due to QCD and the Higgs calculated 
in reference |K| . 

The low energy value of the strong coupling constant considered as an output is given 
by the formula fll4]| , 

a: 1 = a:\M z )\j^ = &z\M z )\ m + ±- - i- In (9) 



*s 



Ait 2tt M z 



The current experimental values of a s at Z-pole extracted from QCD experiments with 
various methods are obtained in Table I || . R refers to the ratio of cross sections or partial 
widths to hadrons versus leptons and the values of strong coupling are in MS renormalization 
scheme. The average value of a s given in Ref. |J is a s (M z ) = 0.118 ± 0.003. 

The values (1),(2) will serve as low energy boundary conditions for the corresponding 
two-loop renormalization group equations. As a high energy boundary conditions for the 
gauge couplings we shall impose unification at a scale Mq, 

di(M G ) = & 2 (M G ) = a 3 (M G ) = a G (10) 

Note that when we depart from the minimal supersymmetric extension of the Standard 
Model and consider SU(5) extensions of it the effect of the superheavy particles with masses 
around the unification scale M G has to be taken into account in the evolution of the gauge 
couplings. 
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TABLE I 


rrocess 


ot s {M z ) 


Deep inelastic scattering 


0.112± 0.006 


R m r lepton decay 


0.122 ± 0.005 


K m I decay 


A 1 AO _l_ A A1 A 

U.lUo ± U.U1U 


Event shapes in e~*~e annihilation 


A 1 094- A AA7 
U. 1ZZ± U.UU 1 


QQ lattice 


0.115± 0.003 


Fragmentation 


0.122± 0.012 


Jets at HERA 


0.121± 0.012 


R in Z° decay (LEP and SLC) 


0.123± 0.006 


Deep inelastic at HERA 


0.120± 0.014 



Table I: Values of a s (Mz) extracted from QCD experiments 



The soft supersymmetry breaking is represented by four parameters M G , Mi/ 2 , A Q and 
B Q of which we shall consider only the first three as input parameters and treat B(Mz), as 
well as the Higgs mixing parameter fi, as determined by the 1-loop minimization equations, 

sin2/3 = - _z BfJ _ 2 (11) 

-(M Z + U ZZ (M Z )) - tan2f3 _ i (12) 
where mf = m\. + /i 2 + — V q^i oov) ■ Note that cases in which the above two equations 

i 

are not satisfied, and thus radiative symmetry breaking does not occur, are rejected. The 
parameter (3 = tan _1 (^) is defined at M Z - M z in (12) denotes the experimental Z-boson 
mass. Thus supersymmetry breaking is parametrized with the input parameters M Q , M\/2, 
A Q , (3{M Z ) and signfi(M z ). We shall take the simplest of boundary conditions at Mq 
assuming universality, 

m,(M G ) = M , M^Mg) = M 1/2 , A t (M G ) = A (13) 



We shall employ the full coupled system of 2-loop renormalization group equations (15 
evolved from M G down to low energies. Since our purpose is to study the effect of high energy 
thresholds in various extensions of SU(5) it will be sufficient to obtain the parameter Ms 
and Ms, appearing in the boundary values of the low energy gauge couplings, by calculating 
the sparticle masses in the step approximation ||16|| . However, we shall introduce finite part 
contributions whenever they are a priori expected to be large |17| , such as QCD corrections 
to the top and gluino masses for instance etc. 

Following the Particle Data Group our basic experimental constraints on supersym- 
metric masses as well as Higgs boson masses are shown in Table II. These limits, impose 
stringent bounds on the extracted values of a s (M z ) and on heavy high energy masses as we 
will see later. 
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TABLE II 

Particle Bound(GeV) 
Neutralinos 

m x o (LSP) >23 

m x ° >52 

m x o >84 

m xl >127 
Charginos 

m x c >45.2 

m x c >99 
Sneutrino 

mp >41.8 
Charged Sleptons 

rrie,jx,T >45 
S quarks 

m 9 - >224 
Gluino 

m-g >154 
Higgs Bosons 

m h >44 
>24.3 

m% >40 
Table II: Experimental bounds on supersymmetric particles and Higgs Bosons 



The purpose of the present paper is to analyze the effect of high energy thresholds on 
unification predictions in the various SU(5) extensions of the MSSM. Our goal is to com- 
plete previous existing analyses which either have not incorporated the constraints imposed 
by radiative symmetry breaking or have not considered the full range allowed for GUT pa- 
rameters. In our consideration of supersymmetric versions of SU(5) we have to take into 
account the constraints imposed by proton decay into K + V IM through dimension 5 operators. 
Assuming that this is the dominant process the proton lifetime is |Tj| [[HJ [^(J 



r ,+ x si , 0.003GW 3 0.67 sin (2/3) M Hc TeV^ 1 ., 

Hp — > K + V u ) = 6.9 x 10 3 V* — 1 — ^ (14) 

(3 n A s (l + y tK ) W^GeV [/(fi, d) + f(u, g)] 

where M Hc is the effective colour triplet mass in GeV,/? n ~ (0.003-0.03)GW 3 , \l+y tK \ > 0.4, 
and A s = (gjgg}) & ( J^g^ ) ~ * ■ The function is defined as 

/(a;,y) = 1 2 (-^ — j" M^r) ~ "T^ — ^M^y)) ( 15 ) 

ar — y z x z — m^j ml- y z — 

and masses are in GeV. The current experimental limit is 

r(p — > K + T7^) > 10 32 yrs (16) 
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In what follows, we choose the most conservative values of (3 n and |1 + y tK \ which are 
0.003GW 3 and 0.4 respectively. 

An additional constraint that will be imposed is the absence of Landau poles on the 
dimensionless couplings of the theory, or equivalently the validity of perturbation theory 
(perturbativity) of those couplings above M G and up to the Planck scale. Although this is 
in general not a severe constraint, it should be taken into account in the cases of extended 
versions of 577(5) due to the existence of a large number of massless particles above M G (see 
for instance Fig. 3c). This constraint is implemented through the numerical solution of the 
1-loop RG equations for the Yukawa couplings of S77(5) above M G (see Appendix). 

2. Minimal SU(5) 



The standard superpotential |2T] in the minimal S77(5) is 



W = -MxTr(S 2 ) + -X^r^ 3 ) + M 2 HH + A 2 77£77 
2 3 

+ V^Y^jH + ^Y^^jH i,j = 1..3 (17) 

S77(5) is spontaneously broken to S77(3) x SU(2) x £7(1) when the adjoint Higgs S develops 
a v.e.v in the direction < £ >= V Diag(2, 2, 2, — 3, — 3). The resulting superheavy masses 
are , 

My = 5g 5 V , M Hc = 5\ 2 V , M s = 5AiV (18) 

We have imposed the usual fine-tuning condition M 2 = 3A 2 V" that gives massless isodoublets 
of 77,77. The mass M-£ stands for the mass of the surviving colour-octet and isotriplet parts 
of E. The leading contribution of these masses to the beta function coefficients of the three 
gauge couplings f = i = 1...3, t = In(^), (63, h, h) = (-3, 1, f ) is 

Ab 3 = 36(Q 2 - M|) + 6(Q 2 - MhJ - A6(Q 2 - My)} (19) 
Ab 2 = 26(Q 2 - Ml) - 66(Q 2 - M v ) (20) 

Ah = h{Q 2 - M 2 H J - 109(Q 2 - M*) (21) 
5 

Demanding perturbativity up to M== ~ 2.4 x 10 18 GeV^ for the couplings appearing in (17), 



/Svr 

we are led after numerically integrating the coupled system of S77(5) renormalization group 
equations, to the inequalities 

Ai < 1.4 ,A 2 < 1.5 ,Y t < 1.5 ,Y h < 1.4 (22) 

at Mq- Note that although in general 

M G = max{M v ,M Hc ,Mx} , (23) 

only the case M G = My can be realized under the combined constraints in our analyses. 
The case in which all superheavy masses are equal although allowed by the bounds given in 
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0.1335 



0.1325 



0.1315 



0.1305 



Minimal SU(5) 

M v =M GUT =max{M v ,M Hi .,M J -) 



\l>0GeV 
A =400 GeV 
' M =750 GeV 
tan$=2.5 
M T =2xlo" GeV 



MSSM 

M w =(8O-180) GeV 



SI 



. 



/ 



s' ^ 



/ 



f 



Proton decav bound 



0.1295 



Minimal SU(5) 

M v =M G1JT =max { M v ,M Hc , MJ 



H"^" Region not allowed by experimental limits 

A =400 GeV J SUSY , ides ' 
1375 «">P=2.J 

M° =3x10" GeV 

i MSXM . — 



2.7 



3.1 



3.5 



3.9 



4.3 



0.1315 - 



0.1305 



0.1295 



M la =(80-170) GeV 
M z =(0.]-2)xl0" GeV 



Region not allowed by proton decay bound 



10 12 



14 



Figure 1: (a),(b) a s (M z ) versus Mq when M Hr , M s , respectively are varied. 



Eq. (22) is the well studied case of the MSSM. In the context of the SU (5) this requires the 
couplings Ai j2 , #5 to be fine tuned according to (18). 

Our standard outputs are the values of the strong coupling, a s (Mz), and the scale Mq 
where the couplings meet. The input values of M , Mu 2 and A Q are always taken to be 
smaller than 1 TeV. It must be noted that in the figures we have chosen input values such 
that the acceptable region in the (a s -Ma) plane is the optimum one in the following sence: 
We vary one parameter at a time while we keep the others constant. This is done for every 
input parameter, until we reach the maximum acceptable area. In the figures shown, we 



have adopted for the mass of the top quark an average value of the CDF and D0 |22 
experimental results, m t = 180 GeV. Note also that, variation of ±5GeV in m t = 180GeV 
results in ±0.0005 and ±0.08 x 10 16 GeV on a s and Mq respectively. In addition, a variation 
of ±0.0007 in the central value of A 7 gives variation of =f0.001 and =f0.5 x 10 16 GeV on a s 
and M G respectively. 

The shaded areas of Figures la and lb represent the allowed parameter space for the 
outputs a s {Mz) and Mq. The results do not depend significantly on M Q and A Q which have 
been chosen to have the representative values shown. The allowed area shrinks with increas- 
ing tan{3 due to the proton decay bound. For smaller values the allowed area shrinks because 
of the perturbativity of Y t (Mo) and the constraint from radiative symmetry breaking. In 
Figure la we have chosen a characteristic value for while we have varied M# c through 
its allowed range of values (1.8 — 4.6) x 10 16 GeV. Analogously in Figure lb we have taken 
Mh c = 3 x 10 16 GeV, while we have varied Ms through the range of values (0.1 — 2) x 10 16 
GeV. The values of a s obtained are in general too large in comparison with the average 
experimental value. Making a parameter search, we conclude that the lowest possible value 
that we can reach in this model is close to ~ 0.130. Nevertheless, there are processes (Table 
I) whose determined a s agrees in a limiting sense with the smallest values in Fig. 1. It 
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should be noted that the effect of high energy thresholds has made the access to the smaller 
values of a s worse than in the case of the MSSM. The general dependence of a s is that it 
increases with increasing M Bc while it decreases with increasing M E . 



3. Missing Doublet Model 

Let us now consider an extended version of the 577(5) known as the Missing Doublet 
Model. This model |23|, constructed in order to avoid the numerical fine-tuning required in 
the minimal SU(5), has instead of the adjoint GUT Higgs a Higgs field in the 75 represen- 
tation as well as an extra pair of Higgses in the 50 + 50 representation. The superpotential 
is 

W = MiTr(£ 2 ) + iAiTr(£ 3 ) + A 2 #£e + A 2 #£e + M 2 ee 

3 

+ ^Y^^+^Y^jH i,j = 1..3 (24) 

Since the 0,6 do not contain any isodoublets, only the coloured triplets obtain masses while 
the isodoublets in H,H stay massless. The v.e.v 

= V{{5 C ) C A {5 C ) D B + 2{5 W ) C A {5 W ) D B - l -5 c A 5 D B - (C <-> D)) (25) 

leads to the masses 

M(3,l,jj) = ^M E ,M(8,1,0) = ^ ,M(8,3,0) = M S 

M{6,2,1) = ,M(1,1,0) = ^M s 

o5 5 

for the remnants of 75. The assignment of the quantum numbers refers to the group SU (3) x 
SU(2) x (7(1) . We shall assume that the parameter M 2 is larger than the GUT scale 
possibly of the order of the Planck mass. Otherwise perturbativity, as can be easily seen, 
cannot be fulfilled. The charge - | colour triplets in H,H and B,B will give one supermassive 
combination of mass M B , — 0(M 2 ) and a light combination of mass 



9 ,32A 2 A 2 wM 



100 v Af ' v M 2 

The mass parameter M E is related to the vector boson mass through 

M E 5Mx (10/3)Ai 1 /5.Ai. 



My 2VTEg 5 V 2VlEg 5 3\ 3 K g 5 



(27) 



The modifications in the beta function coefficients are, 
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Missing Doublet Model 

M V= M GUT= maX i M V. M Hc. 1V y 
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Prolan decay bound 
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4.0 



4.3 
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4.6 



4.9 



Missing Doublet Model 



0.107 



0.106 - 



\l>0 GeV 
A =400 GeV 
M =750 GeV 
ta'n$=2.3 
M„ =3xl0" GeV 



Region nor allowed by exper 'tineiual limits on 
SUSY panicles 



S§ 0.105 - 

, /" *> 
/ 

0.104 



M m =(80-110) GeV 
M z =(O.O5-0.19)xl0"' GeV 



Proton decay bound 



0.103 



Figure 2: (a),(b) a s (M z ) versus Mq when M Hc , M s , respectively are varied. 



A6 3 = -A9(Q 2 - My) + 96>(Q 2 - M|) + 9{Q 2 - 0.8 2 M|) 

+ 10#(Q 2 - 0.4 2 M 2 ) + 36(Q 2 - 0.2 2 M 2 ) + 6(Q 2 - M|J + 6(Q 2 - M^,) 

+ 34fl(Q 2 -M 2 ) (28) 

A6 2 = -6#(Q 2 - M 2 ) + 16fl(Q 2 - M|) + 6#(Q 2 - 0.4 2 M|) 

+ 35fl(Q 2 -M 2 ) (29) 

A6i = -10#(Q 2 -M 2 ) + 10#(Q 2 -0.8 2 M 2 ) + 10#(Q 2 -0.4 2 M 2 ) 

2 ? 1 73 

+ -6(Q 2 - M 2 H J + -6(Q 2 - M 2 Hci ) + — 6(Q 2 -M 2 ) (30) 

Perturbativity above M G , as in the case of minimal SU(5), leads us to the extra constraints 
at M G 

Ai < 0.18 ,Y t < 1.6 ,y 6 < 1.5 (31) 

In this model, we obtain a stronger constraint on Ai, and consequently on M^, due to the 
fact that now we have a larger dimensional representation (75). 

Note that the model as it stands does not contain any /x-term. We assume, however, that 



a /i-term is generated through an independent mechanism [£4[ . 

In Figures 2a and 2b it can be seen that the values of a s obtained are much smaller 
than the average experimental value. We should note that a s is pushed towards smaller 
values due to the splittings within 75 that give a large correction in the opposite direction 
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than in the case of the minimal model [p3 |. The maximum value of a s that we are able to 



obtain in the Missing Doublet Model, is approximately ~ 0.106. As we can see from Table 
I, there are still QCD processes, where the values of a s are in agreement with the results 
of the missing doublet model. The heavy masses M# c and Ms are constrained to be in the 
regions (2.6 - 5.0) x 10 16 GeV and (0.05 - 0.19) x 10 16 GeV respectively. Due to the fact 
that the extracted values of a s in MSSM are greater than 0.125 (for the inputs of Figs 2a,b), 
we do not display in Figures 2a and 2b the corresponding MSSM plane. Experimental limits 
on LSP (see Table II) puts a lower bound on the universal soft gaugino masses such that 
Mi /2 > 80 GeV. 

4. Peccei-Quinn symmetric missing-doublet model 

The problem of proton decay through D=5 operators that is present in the minimal 
SU(5) model provided strong motivation to construct versions of SU(5) with a Peccei-Quinn 
symmetry ||25|| that naturally suppresses these operators by a factor proportional to the ratio 
of the Peccei-Quinn breaking scale over the GUT scale p6| . A Peccei-Quinn version of the 
missing doublet SU(5) model requires the doubling of 5 + 5 and 50 + 50 representations. 
The relevant superpotential terms which must be added to the previous superpotential are, 



\ 2 heo + x 2 hj:q + \' 2 H'ZO' + x' 2 H'm' 

+ M 2 ee' + M' 2 Q'Q + X 3 PH' H' (32) 

P stands for an extra gauge singlet superfield. The charges under the Peccei-Quinn symmetry 
are Of (a/2), <j>({3/2), H(-a), 3T(-2±fi), 6(a), 8(f), &(-a), e'(-^f ), H'(^), if (a), 
P(— (3a + /?)). The breaking of the Peccei-Quinn symmetry can be achieved with a suitable 

gauge singlet system at an intermediate energy < P >= ~ 10 10 — 10 12 GeV. Assuming 
M2, M 2 to be of the order of the Planck scale, we obtain two massive pairs of coloured 
triplets with masses 

M Hc ~ 32A 2 A' 2 ^- , M m ~ 32A 2 A' 2 ^ (33) 

2 2 



somewhat below the GUT scale. Note that My = 2\/l^g^V . In addition we have two pairs 
of isodoublets one of which is massless while the other pair receives the intermediate mass 
Mh' ■ The modifications in the renormalization group beta functions coefficients are, 

Ab 3 = -A6(Q 2 - M v ) + 99(Q 2 - M|) + 6(Q 2 - 0.8 2 M|) 

+ 10#(Q 2 - 0.4 2 M|) + 36(Q 2 - 0.2 2 M|) + 6(Q 2 - M 2 Hc ) + 6(Q 2 - M^J (34) 

Ab 2 = -Q6(Q 2 - M v ) + 1Q6(Q 2 - M|) + Q6(Q 2 - 0.4 2 M|) 

+ 6(Q 2 -M 2 H ,) (35) 

(36) 
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Missing Doublet Model+ Peccei-Quinn 

M v =M GDT ={M v ,M Hc ,M M E } 
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Figure 3: (a),(b),(c) a s {M z ) versus M G when M ff /, M Hc and M s are varied. 



10 



A&i = -10fl(Q 2 - M 2 ) + 10#(Q 2 - 0.8 2 M|) + IW{Q 2 - 0A 2 M 2 ) 

+ \e{Q 2 - M 2 H J + 2 -0{Q 2 - M|J + ^(Q 2 - M 2 H ,) (37) 

Bounds, which come from perturbativity of couplings in this model, are numerically 
similar, to those of the previous one (see Appendix). The extra coupling A3 obeys the 
constraint A3 < 2.7 at Ma- 
lt is evident from figures 3 that the values of a s obtained for the case of this model 
are in excellent agreement with the experiment. The range of these values (0.107 — 0.140), 
covers the experimental average value of a s and lays between the gap of the Minimal SU(5) 
model and the Missing Doublet model. This model possesses an additional parameter, the 
intermediate scale Mjji which increases the values of a s when it takes lower values. The 
allowed range of values for Mi/ 2 has been increased now to 730 GeV or 900 GeV in Figure 
3a and 3b respectively. For this model the allowed range for M can be extended to lower 
values. Figure 3a and 3b have been obtained for M a = 300 GeV. Note however that still 
a s does not depend significantly on M Q . In addition, tan(3 can practically now take much 
larger values, as large as tanf3 ~ 40. Figures 3a, 3b and 3c have been obtained for an 
intermediate value tanj3 = 10. The allowed range of values for the intermediate scale Mh' } 
is (1 - 450) x 10 10 GeV . Similarly 5 x 10 14 < M Hc < M G and 7 x 10 14 < M E < 15 x 10 14 
GeV . Note finally, that the grand unification scale can take values as large as the "string 
scale" for rather large but not excluded values of a s . 



5. A version of SU(5) with light remnants 



Recently there has been some activity around models with gauge group GxG with intent 
to bypass the known problem of k=l superstring constructions where no adjoint Higgs can 
appear in the massless spectrum [[27j. Vector- vector Higgses present in the spectrum of 
GxG can break it into Gdi ag - An SU(5) x SU(5) model with Higgses Z(5,5) + Z(5,5) 
can have renormalizable couplings only of the type ^Z^Z\ = ^Tr(ZZ) to a singlet $, in 
addition to self-couplings of singlets. Thus we could construct an analog SU(5) GUT with 
superpotential, 

W = ^$iTr(S 2 ) + + (38) 

where £ is the adjoint and $1, $2 are singlets. This superpotential is invariant under 
$2 ~~ * — $2- The F-flatness conditions give, apart from < S >= VDiag(2, 2, 2, —3, —3), 

$2 = _^i(30U 2 ) ,$i = (39) 
A2 

With the above superpotential the remnants of E, a colour octet and an isotriplet, stay 
massless. Nevertheless, non-renormalizable terms like 

AW=^Tr(S 4 ) + ^(Tr(S 2 )) 2 (40) 
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Figure 4: (a),(b) a s and M G vs M s in SU(5) model with light remnants. 

can in principle induce a mass of order O(-rr) ~ 10 u GeV or smaller, depending on the 
actual values of A4 and A5. A higher order non-renormalizable term would induce an even 
smaller mass ™ ~ 10 12 GeV. The light S in this model could allow for a large Mq close 
to a string unification scale and thus in such a model there would be no string unification 
mismatch. 

We shall therefore investigate the effects of small Ms on a s and Mq. In order to obtain 
acceptably small values of a s , we choose as input value of Mh c the smallest acceptable one 
as it is shown in Figures 4a and 4b, since an increasing Mh c tends to increase a s . When we 
vary Ms from 10 15 down to 10 13 GeV, we can achieve unification at Mq = 5 x 10 17 GeV 
~ 0(M s t r ing)- Decreasing Ms further towards the intermediate scale, leads to even larger 
values of a s . However, the values of a s are still rather large (> 0.131). This excludes this 
particular version at least in this simple form. 

6. Conclusions 



In this article we have studied various supersymmetric extensions of the Standard Model 
based on the group SU(5). The low energy precision data in conjunction with the existing 
experimental bounds on sparticle masses are known to impose strong constraints if radiative 
breaking of the electroweak symmetry is assumed. There exist several detailed studies in 
the literature in the framework of radiative symmetry breaking |j, [j], |K], TT|, [16|, [IT], which 
however have not considered in detail the effect of the superheavy degrees of freedom included 
in unified schemes. In SUSY GUTs there are additional constraints one has to deal with 
such as the experimental bound on proton's lifetime, the absence of Landau poles beyond 
the unification scale and the appearance of heavy thresholds which influence the evolution 
of the couplings involved. All these affect the low energy predictions. The existing analyses 
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in this direction on the other hand [18, 19, 23, 2C, 26, 27], have not systematically taken 



into account the effect of the low energy thresholds at the level of accuracy required by low 
energy precision data as was done in the previous references. Our analysis combines both 
and takes into account high and low energy thresholds at the accuracy required by precision 
experiments. In particular we have focused our attention on the extracted value of the strong 
coupling constant a s (Mz) , the value of the unification scale Mq, as well as the restrictions 
imposed on the heavy masses in some unifying schemes based on the 577(5). Sample results 
of our findings have been displayed in figures 1-4. 

In the case of the minimal SU(5) we found that the values of the strong coupling constant 
obtained are somewhat larger as compared to the average experimental value of a s (Mz)- 
Also the unification scale Mq differs from the string unification scale by an order of mag- 
nitude if the lower values of a s (Mz) obtained are assumed. Access to small values of the 
strong coupling constant is more difficult than in the MSSM exhibiting the influence of the 
superheavy degrees of freedom in a clear manner. The range of Mh c allowed in this model 
is somewhat limited. The Missing doublet model seems to favour small values of a s (Mz), 
in contrast to those obtained in the MSSM and the minimal SU(5) version. At the same 
time when Mq increases the allowed parameter space shrinks considerably. Proton decay 
along with perturbativity requirements seem to put a stringent constraint on both minimal 
SU(5) and Missing Doublet Model (MDM). In the MDM the large splittings within 75 give 
a high energy threshold effect on a s {Mz) m the opposite direction than in the case of the 
minimal model leading to small values. This is also the case for the Peccei-Quinn version 
of the MDM. However the presence of an extra intermediate scale ameliorates the situation 
allowing to achieve an excellent agreement with the experimental values of a s (Mz)- The 
grand unification scale can take values as large as the string scale at the expense however of 
having rather large values of the strong coupling constant not favoured by all experiments. 
If we consider the range for allowed values of a s (Mz), the minimal SU(5) lays always above 
.130 while MDM lays below .106. The intermediate range between the two can be covered 
by the Peccei-Quinn version of the MDM and coincides with the allowed experimental range. 
Finally, the last considered version, with light S, exhibits unification of couplings at string 
scale but the values of a s obtained are rather large, although within the errors of some 
experiments. 
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Appendix : RGE's above M G 

The Renormalization group equations for the gauge and Yukawa couplings from Mqut 

m p owi o A v i ni8/ 



to ~ 2.4 x 10 18 GeV in the case of minimal SU(5) are pL EU 

Vo7T 11 1 1 1 



3 

dt (4tt 



^ (41) 



^Al ^1 /63 o . 2 ;> 



dt (4tt) 2V 5 



Ai + 3A 2 - 30g 2 5 ) (42) 



dA 2 A 2 ,21 2 53 2 98 2 



A 2 + — A2 - -g£ + 3Y 2 + 4Y 6 ^) (43) 



dt (4tt) 2V 5 1 5 2 5 

dY t Y t ,24 x2 , nv3 , ,,,2 96, 



X 2 2 + 9Y t 2 + AY 2 - -gi) (44) 



dt (4vr) 2V 5 2 * 6 5 

dY » - Y » ,24 x2 . ova , inv2 84 , 



A 2 + 3Y t 2 + 10Y 6 2 - -g 2 ,) (45) 



dt (4tt) 2V 5 2 1 b 5 



where t = ln( — ® t= ] 



The modification of the above system in the case of the missing doublet model is, 

dg 5 17 



dt (4tt) 2 
dAi Ai 56 



9b (46) 
Ai) 2 - 48<fe 2 ) (47) 



dt (4vr) 2V 3 

The other two equations for the top and bottom Yukawa coupling are derived if we set A 2 = 
in the eqs. (43,44) of the minimal model. Due to the fact that we have an extra pair of Higgs 
5 and 5 in the Peccei-Quinn version of the missing doublet model, the only equation that 
changes compared to the missing doublet model is the one for the gauge coupling which 
takes the form 

IF ~ (4^ 5 (48) 

In this model we have an extra coupling A3 whose running is given by the following renor- 
malization group equation 

-dT = jA^ {3 ^-T 9 ^ (49) 
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